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I.  INTRODUCTION 


Guided  waves  in  plates,  which  are  known  as  Lamb  waves,  are  of  considerable  practical  interest 
for  both  nondestructive  evaluation  and  structural  health  monitoring  because  of  their  ability  to 
propagate  relatively  long  distances  while  still  maintaining  sensitivity  to  damage  [1,2].  Because 
of  the  long  propagation  distances,  particularly  as  compared  to  most  bulk  wave  NDE  methods, 
Lamb  wave  are  particularly  sensitive  to  changes  in  the  propagation  environment,  such  as 
temperature,  stress  and  surface  conditions  [3,4,5].  Here  we  consider  the  effect  of  stress  on  Lamb 
waves  propagating  in  homogeneous  and  initially  isotropic  plate  subjected  to  a  homogeneous 
biaxial  stress  field. 

The  theory  to  explain  the  dependence  of  wave  speed  on  stress,  or  acoustoelasticity,  was 
developed  by  Hughes  and  Kelly  by  applying  the  Mumaghan  theory  of  finite  deformation  to 
propagation  of  bulk  elastic  waves  in  an  initially  isotropic  solid  subjected  to  a  static 
predeformation  [6].  They  specifically  consider  the  cases  of  hydrostatic  pressure  and  uniaxial 
compressive  stresses  and  derive  expressions  for  changes  in  shear  and  longitudinal  wave  speeds 
as  a  function  of  applied  stress  for  known  material  properties.  Their  work  was  extended  to 
materials  of  arbitrary  symmetry  by  Toupin  and  Bernstein  [7].  One  motivation  of  this  early  work 
in  acoustoealasticity  was  measurement  of  third  order  elastic  constants;  however,  this  focus  was 
largely  shifted  to  measurement  of  both  applied  and  residual  stresses  [8,9].  The  comprehensive 
article  in  [10]  provides  a  thorough  treatment  of  acoustoelasticity  with  an  emphasis  on  stress 
measurements  via  shear  wave  birefringence. 

Acoustoelastic  Rayleigh  waves  have  also  been  quite  thoroughly  investigated,  e.g.,  [11,12], 
but  there  are  significantly  fewer  published  works  on  acoustoelastic  Lamb  waves.  Husson  [13] 
considered  acoustoelastic  Lamb  waves  from  a  theoretical  point  of  view  and,  among  other  things, 


Acoustoelastic  Lamb  Waves 


Approved  for  public  release;  distribution  unlimited 


Page  2 


predicted  a  strong  frequency  dependence  of  the  acoustoelastic  constants.  Qu  and  Liu  [14] 
generated  dispersion  curves  for  waves  propagating  in  a  stressed  aluminum  plate,  but  did  not 
investigate  dependence  on  the  direction  of  propagation.  Rizzo  and  Lanza  di  Scalea  [15] 
measured  changes  in  wave  speed  with  tensile  loads  for  guided  waves  in  bars,  and  found  a  strong 
frequency  dependence  of  the  acoustoelastic  constants.  Lematre  et  al.  [16]  developed  theory  for 
Lamb  wave  propagation  in  stressed  piezoelectric  plate  structures  but  showed  only  numerical 
results  for  propagation  along  the  direction  of  the  applied  uniaxial  load.  An  in  situ  structural 
health  monitoring  experiment  by  Michaels  et  al.  [17]  clearly  show  changes  in  wave  speed  as  a 
function  of  direction  of  propagation,  but  no  theory  is  presented.  The  semi-analytical  finite 
element  method  has  been  applied  to  acoustoelastic  wave  propagation  in  bar- like  structures  (i.e., 
rails)  [18]. 

In  this  paper  the  theory  of  acoustoelastic  Lamb  wave  propagation  in  a  homogeneous  and 
initially  isotropic  plate  is  developed  for  an  applied  homogeneous  biaxial  stress  by  combining  the 
theories  of  bulk  wave  acoustoelasticity  and  anisotropic  guided  wave  propagation.  The  resulting 
analytical  equations  are  solved  numerically  to  show  the  anisotropic  changes  in  dispersion  for 
various  applied  stresses,  and  these  changes  are  validated  via  experiments  using  an  aluminum 
plate.  Preliminary  results  of  this  present  study  have  been  previously  reported  by  the  authors 
[19,20], 

II.  REVIEW  OF  BULK  WAVE  ACOUSTOELASTICITY 

The  theory  of  bulk  wave  acoustoelasticity  is  reviewed  here  following  the  development  of  Pao 
and  Gamer  [21].  Referring  to  Figure  1,  a  body  is  deformed  from  its  unstressed,  or  natural,  state, 
to  a  statically  deformed,  or  initial,  state;  the  final  state  is  that  of  wave  motion  superposed  on  the 
initial  state.  Coordinates  refer  to  a  material  point  in  the  natural  state,  coordinates  X  refer  to  a 
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material  point  in  the  initial  state,  and  coordinates  x  to  a  material  point  in  the  final  state. 


Deformations  between  the  various  states;  i.e.,  natural  to  initial,  natural  to  final,  and  initial  to 
final,  are  given  by  u1,  uf,  and  u,  respectively: 

u'©  =  X-^, 

uf(£,0  =  x-^,  (1) 

u(§,  t)  =  x  -  X  =  uf  (§,  t)  -  u1  (§). 

To  describe  wave  propagation  in  a  stressed  medium,  the  equation  for  u,  the  incremental 
deformation  between  the  initial  and  final  states  in  a  pre-stressed  medium  must  be  obtained. 

Keeping  the  same  convention  as  is  used  in  [21],  Greek  subscripts  indicate  that  the 
quantities  are  expressed  in  terms  of  the  natural,  or  unstressed  coordinates.  Lagrangian  strain 
tensors  in  the  initial  and  final  states  are, 


E'  =  — 


1  (  du' 


KP  = 


lP  ,  du\  du\ 


2  ^  54  54  54 


l(dua  |  du\  |  du\  du\ 

2  U4  54  54  54 


Note  that  summation  over  repeated  indices  is  implied  here  and  in  subsequent  equations.  If  the 
wave  motion  is  small  compared  to  the  initial  predeformation,  the  incremental  strain  tensor 
between  the  initial  and  final  states  is  given  approximately  by, 


„  <■  if  du„  duB  du)  du ,  du\  du , 

aP  aP  aP  2 1  04  04  04  04  04  04 


If  it  is  assumed  that  the  material  is  hyperelastic,  the  Lagrangian  strain  tensor  can  be  related  to  the 
second  Piola-Kirchoff  stress  tensor  via  a  constitutive  equation;  here  only  the  second  and  third 
order  elastic  constants  are  retained: 
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Tl  —  C  F'  I  1  C  pi  Z7> 

1  ap  ^aPrS^yS  "r  2  ^aPySsn^yd^stj  ’ 

pf  p  pf '  I  1  p  pf  pf 

1  ap  ^aPyS^yS  "r  2  ^aPySsn^yd^stj  • 

Incremental  stresses  and  strains  between  the  initial  and  final  states  are  related  by  subtracting  the 
above  equations  and  discarding  higher  order  terms, 

Tap  ^-'apyS^'yS  ^ ' aPySerf- yS^ eti  •  (^) 

The  infinitesimal  initial  and  incremental  strain  tensors  are, 


-yd 


1 

i 

du  du 

2 

t^6  06] 

en  2 

[06  06] 

(6) 


Equation  (5)  is  thus  the  incremental  constitutive  equation. 

Consider  the  equation  of  equilibrium  for  the  static  predeformation  and  the  equation  of 
motion  for  the  final  state, 


06 


pi  ,  pi  dt/q 

1pa  ^  1py 


06 


=  o, 


06 


rf  +Tf 

1Pa  '  1  py 


<K 

06 


2..f 


=  P 


dzu 


dt 


2  > 


(7) 


where  p  is  the  density  in  the  natural  state.  The  equation  of  motion  for  the  incremental 
displacement  is  obtained  by  subtracting  these  two  equations  and  neglecting  one  higher  order 
term  [21].  The  resulting  equation  of  motion  for  the  incremental  displacement  is, 


dUa 

06 


+  T, 


Pr 


P 


d\ 

dr 


(8) 


Further  simplification  is  possible  by  assuming  that  both  the  material  and  the  static 
predeformation  are  homogeneous,  and  then  substituting  Eq.  (5),  the  incremental  constitutive 
equation,  into  Eq.  (8).  The  resulting  equation  of  motion  for  the  incremental  displacement  is, 
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(9) 


A  =  C  +C  e'  S  +C  _ -  +  C  _ —  +  C  e' 

^afiyS  ^ a Py8  ^ fiSsrj ^ si] u ay  apA8  p^  y~'  Af3y§  p^  aPySsrj^ srj 


Equation  (9)  appears  almost  identical  to  the  usual  elastic  wave  equation,  but  it  should  be  noted 
that  A  does  not  have  the  same  symmetries  as  the  stiffness  tensor  C.  The  only  symmetry  is  Ay^  = 
AkUf,  in  general,  Am  ±  Am  ±  Am. 


III.  THEORY  OF  LAMB  WAVE  ACOUSTOELASTICITY 

Development  of  the  theory  to  explain  the  propagation  of  Lamb  waves  in  pre-stressed  plates 
requires  combining  the  equation  of  motion  and  constitutive  equation  developed  for  bulk  waves  in 
the  previous  section  with  the  theory  for  Lamb  wave  propagation  in  an  anisotropic  plate.  The 
equations  derived  in  the  process  closely  match  those  of  Nayfeh  and  Chimenti  [22]  for  anisotropic 
Lamb  waves. 

The  geometry  and  coordinate  system  used  in  this  paper  are  illustrated  in  Ligure  2. 
Referring  to  the  figure,  the  initial  stresses  are  specified  in  the  primed  coordinate  system,  and  a 
Lamb  wave  propagating  along  any  arbitrary  angle  <p  from  the  xI  direction  is  considered.  The 
unprimed  coordinate  system  is  rotated  through  the  same  angle  <p  to  form  the  primed  coordinate 
system,  and  analysis  is  performed  in  the  unprimed  system  to  simplify  the  mathematics.  Lor  an 
applied  biaxial  stress  along  the  x\  and  x[  directions,  the  initial  stress  tensor  can  be  written  as, 

crn  0  0 

T'=  0  cr22  0  .  (11) 

0  0  0 
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It  is  expressed  in  the  unprimed  system  via  a  rotational  transformation, 


T'  =  B  B  T 

u  I  mu  in  m. 


jn  mn  ’ 


(12) 


where  is  the  cosine  of  the  angle  between  the  x-,  and  the  x'  axes.  For  convenience,  in  this  and 

subsequent  equations,  lower  case  Roman  subscripts  are  used  instead  of  Greek  subscripts  to 
represent  quantities  in  the  natural  coordinate  system.  In  particular,  all  quantities  are  expressed  in 
terms  of  the  natural  coordinates;  for  example,  xj  is  used  instead  of  Q. 

Assuming  that  the  initial  strains  are  small,  the  constitutive  equation  relating  initial  stresses 
and  strains  given  in  Eq.  (4)  can  be  simplified  to, 


T'  =C  e 

1  ij  kl  ■ 


(13) 


For  a  given  applied  stress  field,  the  strains  in  the  unprimed  system  can  be  calculated  by  inverting 
this  equation.  From  this  point  forward,  all  quantities  are  expressed  in  the  unprimed  system. 

The  A  tensor  as  given  by  Eq.  (10)  can  be  simplified  further  by  noting  that  the  displacement 
derivatives  can  be  expressed  in  terms  of  strains  because  for  the  stresses  given  in  Eq.  (11),  the 
rotation  terms  are  zero  for  all  angles  (p.  Thus, 


dulj 

dxk 


(14) 


and 


A...,  =  C..„  +C-,  e  8,  +C  ,e\  +C  ... e 1  +C~W  e  . 

iikl  iikl  limn  mn  ik  iiml  km  mjkl  im  ijklmn  mn  ‘ 


(15) 


The  incremental  stress-strain  relation  of  Eq.  (5)  can  similarly  be  simplified  to 


/  =  B, 


8u,, 


B::1,  —  C;,,,  +  C:-,e,  + 


ij  ijkl  ’  ijkl  ijkl  ^ijml^km  ^ ijklmn*' mn 


(16) 


Acoustoelastic  Lamb  Waves 


Approved  for  public  release;  distribution  unlimited 


Page  7 


Acoustoelastic  Lamb  wave  propagation  for  a  homogeneous  and  biaxial  stress  field  requires 
solving  the  wave  equation  for  the  incremental  displacements  as  given  by  Eq.  (9)  subject  to  stress- 
free  boundary  conditions  at  V3  =  +  d  1 2  with  the  tensor  A  given  by  Eq.  (15)  and  the  stresses  given 
by  Eq.  (16).  This  problem  differs  from  Lamb  wave  propagation  in  anisotropic  media  in  two 
regards:  (1)  as  previously  noted,  the  tensor  A  does  not  have  the  same  symmetries  as  the  stiffness 
tensor  C,  and  (2)  the  boundary  conditions  at  the  plate  surfaces  are  now  written  in  terms  of  the 
incremental  displacements. 

Following  [22],  the  approach  to  solving  this  problem  is  to  assume  solutions  of  the  form, 


_  jj  ^'fOi  +ax3-ct) 


(17) 


where  <f  is  the  wavenumber  in  the  xi  direction,  c  is  the  phase  velocity  along  the  x\  axis  and  a  is 
the  ratio  of  *3  to  x\  wavenumbers.  For  a  specific  value  of  c,  these  solutions  correspond  to 
up-going  and  down-going  bulk  waves  in  the  x\~xi  plane  of  the  plate,  which  are  then  summed 
together  to  form  the  Lamb  wave.  This  approach  is  sometimes  referred  to  as  the  partial  wave 
technique  [23]. 

Substitution  of  Eq.  (17)  into  Eq.  (9)  yields  a  form  of  the  Christoffel  equations, 


Kmn(aWn  =  0, 

where  the  parameters  Kmn  are  given  by, 

K\\=c  p  —  A111  —  a  A313’ 

K 22  —  c~  P  ~  A212  —  CC  A323’ 

K33=c2p°  -Am3-a2A3333, 

K  —  K  —  —  A  —  r/2  A 

1V12  “-21  ^1112  01  ^1323’ 

Kl3  —  K3l  —  —cc{A^i33  +  A^jj), 
K23  —  K32  —  —  a(  Al233  +  A1332). 


(18) 


(19) 
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It  can  be  seen  that  Kmn-Knm.  For  existence  of  non-trivial  solutions  for  the  displacement 
amplitudes  U„,  the  determinant  of  the  K  matrix  must  go  to  zero,  which  produces  a  6th  order 
equation  in  a  with  six  solutions  aq,  q  =  1...6.  The  coefficients  of  the  odd  powers  of  a  are  all 
zero,  resulting  in  a  cubic  equation  in  a  given  by 

P6a6  +  P4a4  +  P2o2  +  P0  =  0,  (20) 


where  the  coefficients  are  listed  in  Appendix  A.  For  a  specific  value  of  the  Lamb  wave  phase 
velocity  c,  solving  this  cubic  equation  yields  six  values  of  a,  which  correspond  to  three  up-going 
bulk  waves  and  three  down-going  bulk  waves. 

The  next  step  is  to  satisfy  the  stress-free  boundary  conditions  on  the  plate  surfaces  as  per 
Eq.  (16).  The  approach  taken  here  is  similar  to  that  taken  in  [22],  and  consists  of  constructing 
displacement  ratios  of  Ui  and  lA  to  U \  for  each  of  the  six  values  of  a. 


U 


V  = 

u, 


2  q 


1? 


u 

w  = 

q  u, 


3  q 


,  q  =  1,2. ..6. 


i  q 


(21) 


Equation  (18)  enables  each  ratio  to  be  expressed  as  a  function  of  the  Kmn  and  the  corresponding 
aq.  The  total  displacement  field  of  the  Lamb  wave  is  the  sum  of  the  displacements  of  the  six 
partial  waves, 


ft 

{u, ,  u2 ,  u3 }  =  ^ {l,  1/  , Wq  }  Ulqe 


i^{xx+aqX2,-ct) 


q=\ 


(22) 


Similarly,  an  expression  for  the  stresses  can  be  derived  by  substituting  Eq.  (22)  into  Eq.  (16). 
The  three  stress  components  in  the  X3  direction  are  then  written  as  a  sum  of  stresses  due  to  the  six 
individual  bulk  waves, 


4(xi+a„x3-ct) 


(23) 
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where  the  parameters  Dmq  are  given  by, 


DXq  =  B33n  +  B33l2Vq  +  ccqB3333Wq, 

E2q  —  aq  (-®1313  +  E\323^q  )  +  ^1331^9  ’  (24) 

E3 q  ~  aq  ( ^1 323  +  ^2323 Vq  )  +  ^1332^? * 

Applying  the  stress-free  boundary  conditions  at  the  surface  of  the  plate  requires  setting  T 13,  T2 3 
and  Z33  to  zero  at  x2  =  ±d/ 2,  and  yields  six  equations  in  terms  of  the  six  displacement 
amplitudes,  U \q,  of  the  six  partial  waves,  q  =  1,2... 6.  The  determinant  of  coefficients  must  go  to 
zero  to  obtain  nontrivial  solutions  for  these  six  displacement  amplitudes.  This  determinant  is, 


DnEi 

d12e2 

DnE3 

D\4E4 

Dl5E5 

E16E6 

E*1\E\ 

D22E2 

D23E 3 

E24E4 

E25E5 

E26E6 

D3A 

D32E2 

D33E3 

E34E4 

E35E5 

E36E6 

=  o, 

DnEi 

D\2E2 

D13E3 

E14E4 

D\A 

E16E6 

E^1\E\ 

D22E2 

D23E3 

E24E4 

E25E5 

E26E6 

D31Ex 

D32E2 

D33E3 

E34E4 

E35E5 

E36E6 

where  E  =  e‘,ciqdn  and  E  =e  After  several  number  of  row  and  column  operations, 

Eq.  (25)  decouples  into  two  equations, 

fs  ( co ,  c)  =DnGx  co  t(yax )  +  D13G3  cot  (ya3 )  +  DX5G5  cot(ya5 )  =  0,  (26) 

fa  (ox  c )  =  DUGX  tan(;/o'l )  +  D13G3  tan(r«3 )  +  DI5G5  tan(ya5 )  =  0.  (27) 

The  first  equation  corresponds  to  the  symmetric  modes  and  the  second  one  corresponds  to  the 
antisymmetric  modes.  In  both  equations  y  =  qd/2  =  cod/ (2c),  where  co  is  the  angular  frequency. 
The  parameters  G,„  are  given  as 

Gl  =  D23D35  ~  D33D25  ’  G3  =  D3lD25  ~  D2\D35  ’  GS  =  D2lD33  ~  D3\D23  •  (28) 
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Solving  Eqs.  (26)  and  (27)  yields  die  dispersion  curves  relating  phase  velocity  and  angular 
frequency  for  the  symmetric  and  antisymmetric  Lamb  wave  modes,  respectively.  It  should  be 
noted  that  solving  these  equations  is  not  trivial;  a  brief  description  of  the  numerical  method  is 
given  in  Appendix  B. 

IV.  SELECTED  ANALYTICAL  RESULTS 

This  section  presents  dispersion  curves  generated  by  solving  Eqs.  (26)  and  (27)  numerically  for 
specific  materials.  The  first  material  considered  is  6061-T6  aluminum,  which  was  selected  to 
match  subsequent  experiments.  Material  constants  used  are  listed  in  Table  1,  and  the  third  order 
elastic  constants  were  taken  from  [24].  It  should  be  noted  that  although  third  order  elastic 
constants  are  available  in  the  literature  for  several  aluminum  alloys  [24,25,26],  their  values  are 
quite  variable.  It  is  not  clear  whether  the  variability  is  caused  by  differences  in  experimental 
methods  or  actual  material  variations.  The  second  material  considered  was  steel  because  of  its 
widespread  use.  Unfortunately  third  order  elastic  constants  for  steel  are  not  widely  reported;  the 
ones  used  here  are  summarized  in  Table  2  and  are  for  Hecla  37  steel  [26]. 

Figures  3(a)  and  3(b)  present  the  family  of  symmetric  and  antisymmetric  dispersion  curves 
for  aluminum  under  a  uniaxial  stress  of  an  =  100  MPa  and  for  waves  propagating  at  an  angle  of 
45°  to  the  stress  direction;  note  that  the  SHo  mode  is  not  shown.  To  illustrate  the  variation  of  the 
function  values  as  they  go  to  zero,  log  I  fs  (a>,  c)  I  and  log  I  fa(co,c)  I  have  been  plotted  in  the 
background  for  both  plots.  These  dispersion  curves  appear  virtually  identical  to  the  ones 
obtained  using  the  Rayleigh-Lamb  equations  for  the  no-load  case  since  the  velocity  changes  due 
to  the  applied  stress  are  not  noticeable  at  this  scale.  These  and  all  subsequent  curves  refer  to 
phase  velocities  computed  in  the  natural,  or  undeformed,  coordinate  system,  which  is  consistent 
with  the  derivations  of  Section  III. 
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Figure  4  presents  the  variation  of  phase  velocity  of  the  So  mode  at  250  kHz  with  angle  for  a 
range  of  biaxial  stresses  and  for  a  plate  thickness  of  6.35  mm.  Note  that  the  load  is  applied  along 
the  y  axis  so  that  waves  propagating  along  the  direction  of  applied  load  are  at  90°.  An  interesting 
prediction  made  by  this  plot  is  that  the  angle  at  which  the  phase  velocity  change  is  zero  (about 
27°)  is  independent  of  stress.  It  can  also  be  seen  that  at  a  specific  angle,  the  change  of  phase 
velocity  is  linear  with  load. 

Figure  5  compares  phase  velocity  changes  as  a  function  of  frequency  for  the  So  mode  at 
varying  unixial  stresses  an  for  waves  propagating  at  an  angle  of  45°.  This  plot  shows  large 
changes  in  phase  velocity  for  lower  frequencies  that  reduce  and  flatten  in  the  high  frequency 
region.  It  is  evident  that  the  change  of  phase  velocity  is  linear  with  load  at  a  specific  frequency. 

Figures  6(a)  through  6(d)  compare  the  phase  velocity  change  with  frequency  at  different 
angles  for  the  So,  Si,  Ao  and  A|  Lamb  wave  modes,  respectively.  The  plot  for  A0  is  particular 
interesting  in  that  it  predicts  isotropic  phase  velocities  for  the  Ao  mode  for  the  frequency- 
thickness  product  of  approximately  187  MHz-mm,  although  at  this  frequency  the  velocity 
change  from  the  unstressed  condition  is  non-zero.  This  phenomenon  should  not  be  confused  as 
isotropic  wave  propagation  since  the  group  velocity  at  this  frequency  is  certainly  angle 
dependent. 

Figure  7  is  a  plot  of  the  So  and  SHo  modes  under  the  influence  of  applied  uniaxial  stress 
and  shows  interesting  behavior.  The  So  and  SHo  modes  split  and  form  two  new  continuous 
modes  at  their  original  point  of  intersection  (from  the  unstressed  case)  when  a  stress  is  applied. 
Because  of  symmetry,  this  splitting  does  not  occur  at  either  0°  or  90°,  and  is  similar  to  that  which 
also  occurs  for  anisotropic  materials  [27,28].  In  the  earlier  plots,  this  splitting  was  not  shown  to 
facilitate  comparison  of  the  perturbed  mode  to  the  So  mode  of  the  unstressed  material. 
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Figure  8  shows  the  Ao  mode  for  steel  under  a  uniaxial  stress  of  an  =200  MPa.  Steel  also 
shows  the  isotropic  phase  velocities  at  lower  frequencies  for  the  Ao  mode  at  a  frequency- 
thickness  product  of  approximately  0.35  MHz-mm.  Additional  results  for  steel  are  not  shown 
because  they  are  qualitatively  the  same  as  for  aluminum. 

V.  EXPERIMENTAL  RESULTS 

An  experiment  was  performed  to  measure  changes  in  phase  velocity  as  a  function  of  applied  load 
for  various  Lamb  modes  propagating  in  an  aluminum  plate  with  dimensions  of  305  mm  x  610 
mm  x  6.35  mm.  The  plate  was  instrumented  with  ten  PZT  discs  as  shown  in  Figure  9,  and  the 
transmitting  transducers  were  excited  by  a  five-cycle  Hanning  windowed  sinusoid.  The 
following  three  frequencies  were  utilized  to  produce  an  unambiguous  first  arrival  of  the  indicated 
Lamb  wave  mode:  250  kHz  (So  mode),  400  kHz  (Ai  mode),  and  600  kHz  (Si  mode).  Signals 
were  recorded  from  9  of  the  45  possible  transmit-receive  pairs  as  shown  in  Figure  9(b)  at  eleven 
uniaxial  loads  (0  MPa  to  57.5  MPa  in  steps  of  5.75  MPa). 

Time  shifts  were  measured  by  identifying  a  zero  crossing  in  the  center  of  each  direct  arrival 
echoes,  and  then  tracking  that  zero  crossing  as  a  function  of  applied  uniaxial  load.  For  a 
particular  time  shift,  the  corresponding  change  in  phase  velocity  can  be  calculated  as  [19], 

c2At 

Ac'=-i r-  (29) 

Although  there  is  also  a  change  of  distance,  it  does  not  appear  in  Eq.  (28)  because  all 
computations  are  in  the  natural,  or  undeformed,  coordinate  system. 

Figure  10(a)  shows  the  linear  variation  of  the  changes  in  phase  velocity  with  stress  for 
propagation  of  the  three  Lamb  wave  modes  along  the  direction  of  applied  stress  (the  y  axis).  The 
results  are  further  corroborated  by  experimental  measurements  using  transducer  #1  as  a 
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transmitter  and  transducer  #6  as  a  receiver,  which  show  similar  linear  changes  in  phase  velocity. 
Figure  10(b)  shows  the  variation  of  phase  velocity  with  angle  for  the  case  of  uniaxial  stress  of 
a22  =57.5  MPa.  The  experimental  values  were  obtained  by  fitting  a  line  to  the  time  shift  versus 
load  curve,  and  then  calculating  the  time  shift  at  57.5  MPa  from  the  linear  fit.  Both  theoretical 
and  experimental  data  closely  follow  a  sinusoidal  profile.  The  systematic  differences  between 
theory  and  experiment  are  probably  due  to  the  well-known  difficulties  in  accurately  obtaining 
third  order  elastic  constants.  The  linear  dependence  of  phase  velocity  changes  with  load  and 
sinusoidal  dependence  with  angle  of  propagation  is  expected  since  that  is  the  case  for  bulk 
waves,  but  it  is  reassuring  that  both  characteristics  are  validated  by  both  numerical  and 
experimental  results. 

Figure  11,  which  us  analogous  to  Figure  12(b),  shows  the  variation  of  phase  velocity  with 
angle  but  with  the  theoretical  curves  calculated  using  different  (but  plausible)  values  for  the 
TOECs.  For  these  values,  selected  as  l  =  -181  MPa,  m  =  -289  MPa,  and  n  =  -336  MPa,  the  So 
and  Ai  theoretical  curves  are  in  excellent  agreement  with  the  experimental  data,  and  the  Si  curve 
is  in  reasonable  agreement.  This  significant  improvement  not  only  supports  the  speculation  that 
the  deviation  between  theory  and  experiment  is  largely  explained  by  errors  in  TOECs,  it  also 
suggests  that  guided  wave  measurements  could  be  used  to  estimate  TOECs. 

VI.  SUMMARY  AND  CONCLUSIONS 

Described  here  is  the  theory  for  acoustoelastic  Lamb  wave  propagation,  and  both  numerical  and 
experimental  results  are  shown  for  the  multiple  Lamb  wave  modes.  Unlike  bulk  waves,  the 
acoustoelastic  behavior  is  frequency  dependent.  Experimental  results  that  measure  changes  in 
phase  velocity  as  a  function  of  both  load  and  propagation  angle  agree  reasonably  well  with 
theory.  Both  numerical  and  experimental  data  show  the  expected  linear  dependence  of  phase 
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velocity  changes  with  load  and  sinusoidal  dependence  with  angle  of  propagation,  all  at  a  single 


frequency.  However,  systematic  differences  between  theory  and  experiment  clearly  indicate  that 
additional  measurements  should  be  performed  to  independently  measure  third  order  elastic 
constants  for  the  specific  material  under  consideration. 

The  robustness  of  the  Lamb  wave  results  as  a  function  of  propagation  angle  suggest  that 
Lamb  wave  measurements  may  provide  an  alternative  means  of  measuring  third  order  elastic 
constants.  This  would  be  challenging  because  unlike  bulk  wave  measurements,  there  is  no 
analytical  expression  relating  third  order  elastic  constants  to  measured  time  shifts,  much  less  a 
closed  form  solution  for  the  inverse  problem.  Thus,  any  means  of  using  Lamb  wave  data  to 
obtain  third  order  elastic  constants  would  have  to  be  numerical  in  nature  with  the  resulting  issues 
of  existence  and  uniqueness. 

Additional  experiments  are  also  needed  to  verify  some  of  the  interesting  theoretical 
predictions,  particularly  those  shown  in  Figure  4(c)  regarding  isotropic  phase  velocities  changes 
in  the  low  frequency  region  for  the  Ao  mode.  Future  theoretical  work  should  consider  extension 
of  this  theory  to  materials  of  higher  symmetry,  which  is  believed  to  be  straightforward  but 
tedious.  In  terms  of  numerical  work,  there  are  certain  regions  in  the  anti- symmetric  spectrum 
that  may  indicate  the  presence  of  roots  (SHo  like  solutions),  although  they  do  not  correspond  to 
any  mode  for  the  unstressed  case.  Most  likely  these  are  false  roots  are  caused  by  imperfections 
in  the  numerical  method,  but  they  need  further  investigation. 
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APPENDIX  A.  Polynomial  Coefficients 

Coefficients  of  the  6th  order  polynomial  in  Eq.  (20)  are  as  follows: 


P(,  (A323  A3nA323)A333 

p  -  A2  A  +  A  A2  _9  A  A  A  _  O  A  A  A  +  4  A2 

1  4  ~  ^1233^1313  T  ^13 13^1323  Z,/11 133^1323^1332  z"/11323'rl1331'rl1332  T  '211313'rl1332 

+  2A1233(— A1323(A1133  +  A1331)  +  A1313A1332)  +  A1133A2323  — A1313A2323 
+  2A1133A1331A2323  +  Ami^2223  _  (A212A313  _  ^All2A323  +  A  1 11 A323 )A 
+  c  Ai(— A323  A313A323  (A313  A323  )■ A333 ) 

p  -  A2  A  +  4  A2  _  A  A2  +94  4  4  -94  4  4 

J2  ^1133^1212  T  ^1 111^*1233  'rl1212'rl1313  T  z"rl1112'rl1313'rl1323  z"rl1112'rl1233'rl1331 

+  4  4 2  +94  4  4  -94  4  4  +4  42 

T ^1212^1331  T  ^^1 1 1 1^1233^1332  ^^l  1 12^1331^1332  T  ^1 1 1 1^1332 

— 2Au33(— Ai212A133i  +Au12(A1233  +  A1332))  —  AnilA1313A2323 

+(A  112  A 1 1 1 A212  ) Al333  C  A)  (A313  A323  A333) 

+C  Pq(-~ 2A112A323  (A  133  A331  )  _  (A233  A332) 

A313  ^  A212  +  A313  +  A323  )  "*"  A212  A333  Alll(A>323  A333)) 

Pq  — (— A112  (c  Po  —  Alll)(C  Po  ~  A212))(C  Po  —  A313) 


3333 


(Al) 


Appendix  B.  Numerical  Method 

This  appendix  summarizes  the  steps  of  the  numerical  method  for  solving  Eqs.  (26)  and  (27)  for 
the  symmetric  and  antisymmetric  modes. 


The  first  step  is  to  pick  a  plate  wave  velocity  c  and  find  the  corresponding  aq  from  the 
polynomial  in  Eq.  (20).  Next,  express  Eq.  (25)  as 

fs  ( co ,  c)  =  Dl  j  G,  (a3 ,  a- )  cot(y«, )  +  D,,G3  (a, ,  a5 )  cot  (ya3 )  +  Dl5G5  (a{ ,  a, )  cot(ya5 ) 

(B 1) 

=  Dx  1Hl  (a, ,  a3  ,a5')  +  DUH3  (cq ,  a3 ,  a5 )  +  Dl5H5  (a, ,  a3 ,  a5 )  =  0. 

For  real  aq,  all  Ds,  Gs  and  cotangent  functions  are  real  and  therefore /s  is  real.  Diq  is  even  in  aq, 
Di q  and  D3q  are  odd  in  aq,  which  implies  Gn(ap,  aq)  is  odd  in  both  ap  and  oq.  Since  cot(yaq)  is 
odd  in  aq,  Hn(ap,  aq,  ar)  is  odd  in  all  three  a.  For  odd  functions,  the  Taylor  series  has  only  odd 
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terms  and  this  implies  that  for  imaginary  arguments,  the  value  of  the  odd  functions  is  always 
imaginary.  By  carefully  analyzing  the  symmetries  for  all  terms  it  can  be  seen  that  for  a  mixture 
of  real  and  imaginary  as,/s  has  a  value  that  is  purely  imaginary  or  purely  real. 

Next,  based  on  the  values  of  a,  we  pick  one  of  two  root  finding  algorithms: 

(1)  For  values  of  c  that  produce  purely  imaginary  and  real  values  of  a: 

•  Sweep  across  co  and  find  pairs  of  co  values  where  the  value  of/s  changes  sign  (this  is  not  a 
problem  since /s  values  at  all  co  for  this  value  of  c  are  real  or  all  imaginary).  Each  of  the 
pairs  thus  obtained  mark  intervals  in  which  solutions  to  the  equation  fs  =  0  are  guaranteed 
to  exist  (and  these  intervals  can  be  arbitrary  small).  These  (co,  c)  pairs  constitute  points 
on  the  dispersion  curves. 

•  Any  real  value  of  a  will  produce  infinities  in  the  cotangent  functions.  Since  a  sign 
change  occurs  at  aninfinity,  these  are  actually  false  roots  that  must  be  handled  by 
excluding  pairs  of  co  values  that  correspond  to  multiples  of  rnt  in  the  cotangent  functions. 
Cotangents  with  imaginary  values  of  a  do  not  produce  any  infinities. 

(2)  For  values  of  c  with  any  complex  values  of  a: 

•  Ds  and  Gs  are  complex  valued  and  therefore  fs  is  complex  valued.  Sweep  across  co  for 
pairs  of  co  where  both  the  real  and  imaginary  components  of fs  switch  sign  together.  Also 
project  real  and  imaginary  components  of  the  function  value  fs  on  the  6  =  45°  and  9  = 
135°  lines  and  check  for  sign  changes  in  both  projected  directions.  This  procedure  finds 
roots  that  occur  when  the  curve  of  /s  as  a  function  of  co  is  tangential  to  either  the  real  or 
imaginary  axis.  Either  a  sign  change  in  both  real  and  imaginary  components  or  ones  in 
both  projected  values  yields  an  interval  with  a  solution.  In  general  for  this  case,  it  is 
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possible  to  get  false  solutions  but  these  can  be  discarded  later.  Cotangents  corresponding 
to  complex  values  of  a  do  not  produce  any  infinities. 

Finally,  points  that  are  nearby  in  the  set  of  solutions  are  connected  to  obtain  dispersion  curves. 
Ambiguous  regions  can  be  scanned  with  higher  resolutions  to  obtain  more  points.  The  exact 
same  method  applies  to  the  function  for  antisymmetric  modes /a. 
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Table  1.  Material  constants  for  6061-T6  aluminum  [24]. 


Parameter 

Value 

Units 

2 

54.308 

GPa 

P 

27.174 

GPa 

1 

-281.5 

GPa 

m 

-339.0 

GPa 

n 

-416.0 

GPa 

Po 

2704 

kg/m3 

Table  2.  Material  constants  for  Hecla  37  steel  [26]. 


Parameter 

Value 

Units 

2 

111 

GPa 

P 

821 

GPa 

1 

-461 

GPa 

m 

-636 

GPa 

n 

-708 

GPa 

Po 

7823 

kg/m3 
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Figure  1.  Deformation  of  a  body  from  its  natural  (undeformed)  state  \  to  an  initial  state  of  static 
deformation  X  to  a  final  state  of  wave  motion  x. 

Figure  2.  Geometry  for  Lamb  wave  propagation  in  a  pre-stressed  plate.  Stresses  are  applied 
along  the  principal  directions  in  the  primed  coordinate  system,  and  Lamb  waves  propagate  along 
the  x\  axis. 

Figure  3.  (color  online)  Dispersion  curves  for  waves  propagating  in  aluminum  at  an  angle  of  45° 
to  an  applied  uniaxial  applied  stress  of  100  MPa.  (a)  Symmetric  modes,  and  (b)  anti- symmetric 
modes. 

Figure  4.  (color  online)  Change  in  phase  velocity  versus  propagation  angle  for  the  So  mode  in 
aluminum  (thickness  of  6.35  mm)  at  250  kHz  and  for  different  uniaxial  loads  applied  along  the  y 
axis  (90°). 

Figure  5.  (color  online)  Change  in  phase  velocity  versus  frequency  for  the  So  mode  propagating 
in  aluminum  at  an  angle  of  45°  to  various  uniaxial  loads  applied  along  the  x  axis. 

Figure  6.  (color  online)  Changes  in  phase  velocity  versus  frequency  for  waves  propagation  in 
aluminum  at  various  angles  to  a  100  MPa  uniaxial  load  applied  at  0°.  (a)  So,  (b)  Si,  (c)  Ao,  and 
(d)  Ai  mode. 

Figure  7.  (color  online)  Dispersion  curves  for  the  So  and  SHo  modes  for  a  uniaxial  load  of  100 
MPa  and  propagation  at  an  angle  of  45°  to  the  applied  load.  The  inset  plot  shows  that  the  modes 
are  mixed  and  do  not  cross. 
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Figure  8.  (color  online)  Changes  in  phase  velocity  for  the  Ao  mode  propagating  in  steel  for  a 
uniaxial  load  of  100  MPa  applied  along  the  x  axis  (0°).  The  inset  plot  shows  that  the  changes  are 
isotropic  for  a  frequency-thickness  product  of  approximately  0.35  MHz-mm. 

Figure  9.  (a)  Photograph  and  (b)  sketch  of  transducers  mounted  to  a  6.35  mm  thick  aluminum 
plate.  The  sketch  shows  the  nine  propagation  paths  used  to  characterize  acoustoelastic  Lamb 
wave  propagation. 

Figure  10.  (color  online)  Comparison  of  theory  and  experiment  for  propagation  in  aluminum 
(thickness  of  6.35  mm)  with  an  applied  uniaxial  load  at  90°.  (a)  Changes  in  phase  velocity 
versus  load  for  waves  propagating  at  an  angle  of  90°,  and  (b)  changes  in  phase  velocity  versus 
propagation  angle  for  an  applied  load  of  57.5  MPa. 

Figure  11.  (color  online)  Comparison  of  theory  and  experiment  with  an  applied  uniaxial  load  of 
57.5  GPa  at  90°  using  modified  TOECs  (/  =  -181  MPa,  m  =  -289  MPa,  and  n  =  -336  MPa). 
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Figure  1.  Deformation  of  a  body  from  its  natural  (undeformed)  state  \  to  an  initial  state  of  static 
deformation  X  to  a  final  state  of  wave  motion  x. 


Figure  2.  Geometry  for  Lamb  wave  propagation  in  a  pre-stressed  plate.  Stresses  are  applied  along  the 
principal  directions  in  the  primed  coordinate  system,  and  Lamb  waves  propagate  along  the  xi  axis. 


Acoustoelastic  Lamb  Waves 


Page  25 


Approved  for  public  release;  distribution  unlimited 


(a) 


(b) 


Figure  3.  (color  online)  Dispersion  curves  for  waves  propagating  in  aluminum  at  an  angle  of  45°  to  an 
applied  uniaxial  applied  stress  of  100  MPa.  (a)  Symmetric  modes,  and  (b)  anti-symmetric  modes. 
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Figure  4.  (color  online)  Change  in  phase  velocity  versus  propagation  angle  for  the  So  mode  in  aluminum 
(thickness  of  6.35  mm)  at  250  kHz  and  for  different  uniaxial  loads  applied  along  the  y  axis  (90°). 
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Figure  5.  (color  online)  Change  in  phase  velocity  versus  frequency  for  the  S0  mode  propagating  in 
aluminum  at  an  angle  of  45°  to  various  uniaxial  loads  applied  along  the  x  axis. 
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Figure  6.  (color  online)  Changes  in  phase  velocity  versus  frequency  for  waves  propagation  in  aluminum 
at  various  angles  to  a  100  MPa  uniaxial  load  applied  at  0°.  (a)  S0,  (b)  Si,  (c)  A0,  and  (d)  A,  mode. 
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Figure  7.  (color  online)  Dispersion  curves  for  the  S0  and  SH0  modes  for  a  uniaxial  load  of  100  MPa  and 
propagation  at  an  angle  of  45°  to  the  applied  load.  The  inset  plot  shows  that  the  modes  are  mixed  and  do 
not  cross. 
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Figure  8.  (color  online)  Changes  in  phase  velocity  for  the  A0  mode  propagating  in  steel  for  a  uniaxial 
load  of  100  MPa  applied  along  the  x  axis  (0°).  The  inset  plot  shows  that  the  changes  are  isotropic  for  a 
frequency-thickness  product  of  approximately  0.35  MHz-mm. 
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Figure  9.  (a)  Photograph  and  (b)  sketch  of  transducers  mounted  to  a  6.35  mm  thick  aluminum  plate.  The 
sketch  shows  the  nine  propagation  paths  used  to  experimentally  characterize  acoustoelastic  Lamb  wave 
propagation. 
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Figure  10.  (color  online)  Comparison  of  theory  and  experiment  for  propagation  in  aluminum  (thickness 
of  6.35  mm)  with  an  applied  uniaxial  load  at  90°.  (a)  Changes  in  phase  velocity  versus  load  for  waves 
propagating  at  an  angle  of  90°,  and  (b)  changes  in  phase  velocity  versus  propagation  angle  for  an  applied 
load  of  57.5  MPa. 
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Figure  11.  (color  online)  Comparison  of  theory  and  experiment  with  an  applied  uniaxial  load  of  57.5  GPa 
at  90°  using  modified  TOECs  (/  =  -181  MPa,  m  =  -289  MPa,  and  n  =  -336  MPa). 
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